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We develop a practical Hartree-Fock theory for trapped Bose and Fermi gases that interact with dipole-dipole 
interactions. This theory is applicable at zero and finite temperature. Our approach is based on the introduction 
of local momentum distortion fields that characterize the exchange effects in terms of a local effective potential. 
We validate our theory against existing theories, finding excellent agreement with full Hartree-Fock calculations. 
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Introduction: Phenomenal progress in the production of 
ultra-cold quantum gases with magnetic 0J0I and electric 
|H]|6l dipoles has opened up an important new manybody sys- 
tem 0. The key feature of these gases is that the constituent 
particles interact via a dipole-dipole interaction (DDI) that is 
long-ranged and anisotropic. 

There has been considerable success in the development of 
theory for dipolar Bose-Einstein condensates, in which all the 
atoms occupy a single mode that is described by the mean- 
field Gross-Pitaevskii equation 0. However in situations 
where many modes are occupied (i.e. a Bose gas at finite tem- 
perature or a Fermi gas) the meanfield treatment of the non- 
local exchange interaction is technically challenging. This is- 
sue is most pronounced in the experimentally relevant case of 
trapped samples where both direct and exchange effects con- 
tribute. To date, calculations including exchange have been 
performed by two groups for normal Bose and Fermi gases 
within the Hartree-Fock (HF) approximation l9 ffl"2ll . and for 
small quasi-two dimensional condensates lfl3ll within the HF- 
Bogoliubov-Popov approach. These calculations are numer- 
ically intensive and are only practical in cases where the di- 
mensionality is reduced, either through cylindrical symmetry 
or by tight confinement. Some simple variational lfl4l and 
phenomenological IfTD treatments of exchange have been in- 
vestigated (see comparisons to HF calculations in [9—12 1). We 
also note the application of beyond-meanfield Monte-Carlo 
methods to two-dimensional gases fl6l . 

Exchange effects are predicted to cause dipolar gases to un- 
dergo momentum space magnetostriction [12], and have a sig- 
nificant role in mechanical stability |9, 10 1. A number of stud- 
ies of homogeneous Fermi systems have shown the impor- 
tance of exchange for various phase transitions and other phe- 
nomena [17 1, but extending these predictions to the trapped 
system remains an outstanding problem. 

Here we report on the development of a tractable Hartree 
Local-Fock (HLF) theory for trapped dipolar gases that ac- 
curately describes both direct and exchange interactions. Our 
theory is based on the semiclassical HF approximation (avoid- 
ing the need to diagonalize for modes), a theory that has been 
extensively applied to gases with contact interactions |fl~8), 
and provides a good description of experiments (e.g. see lfT9l ). 
The HLF theory is derived by introducing a pair momentum 
distortion fields that simplify the exchange term to a local po- 
tential. This approach provides insight into the manifestation 
of exchange interactions in the dipolar gas, and opens a path 
for developing meanfield theories in the superfiuid regime. 



We validate the HLF theory against HF and Hartree cal- 
culations for Bose and Fermi systems at zero and finite tem- 
perature. The HLF theory is vastly faster and more resource 
efficient: a HF calculation taking 40 hours is reduced to 2 sec- 
onds with HLF0 

System: We consider a gas of spin polarized particles that 
interact by a DDI of the form 



tfdd(r) 



C dd 1 - 3cos 2 6> 
17 IrP ' 



(1) 



where Cdd = MoMm f° r magnetic dipoles of strength /j, m and 
d 2 /eo for electric dipoles of strength d, and is the angle be- 
tween the dipole separation r and the polarization axis, which 
we take to be the z direction. The particles also interact via a 
contact interaction of strength g (note g = for spin-polarized 
fermions) and are taken to be confined within a trap I7t r (x) of 
arbitrary geometry. 

HLF theory: The single particle Wigner distribution func- 
tion, within the semiclassical approximation, is given by 



W(x,k) 



1 



,P[e(x,k)-ii] 



(2) 



where r\ — 1 for bosons and r/ = — 1 for fermions, /i is the 
chemical potential, and (3 = l/fe^T is the inverse tempera- 
ture. The HLF theory is based on a trial dispersion relation 



e(x, k) = — [K p (x)k; + k z (x)K\ + V e s(x), (3) 



where k„ 



k 2 



v p — y n x t ky, and the effective potential V q (x), 

as we show below, includes the influence of trap, direct and 
exchange interaction. We have also introduced local momen- 
tum distortion fields k p (x) and k z (x), which describe a spa- 
tially varying anisotropy of the momentum distribution with 
respect to the z axis (i.e. direction of dipole polarization). Be- 
cause the momentum distortion determines the anisotropy of 
the pair correlation function [ 12 1, these fields parameterize the 



1 Assuming a cylindrically symmetric trap to make HF calculations feasible, 
with (N P ,N Z ) grid points in (radial, axial) directions (both momentum 
and spatial directions for HF and just spatial directions for HLF), for HF 
calculations the slow step is calculating <I>£;(x, k) which is 0[(N p N z ) 3 ], 
and for HLF the slow step is calculating <5(x) which is 0(N p N z ). The 



example times given are for fixed fj, with N p 



N, 



100. 
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exchange interaction within the HLF theory. We note that the 
cylindrical symmetry of the DDI ([TJ allows us to make the de- 
composition into K p (x) and K z (x) fields, irrespective of the 
trap geometry. 

Using the trial dispersion the position density is given by 



/dV 
(2^ (x ' k) 



^b(x) 



(4) 



where AdB ( x ) = \j2-n:[3h 2 /m* (x) is the thermal de Broglie 
wavelength with spatially dependent effective mass m* (x) = 

m/ [ Kp (x) 2 Kz (x)] 1/3 , and Q{z) = £~ x ^"V/*" is the 
polylogarithm function. 

By applying a variational principle to the free energy, we 
derive the equations for the local momentum distortion fields 
and the effective potential which define the HLF theory. The 
exact equilibrium (grand) free energy f2 cx satisfies ll20l 

Hex < fW = °o - E + (H) , (5) 
where ^hlf is the HLF free energy and 

Eo = j ^e(x,k)W(x,k), 



(6) 
(7) 



= /dx[Jf(x)+V rf (x)n(x)] 



(8) 



are the free energy and single particle energy, respectively, 
with 



X(x) = 



3 kBT /-n f P P[»-VMx)} 
2A dB (x)W 



(9) 



The quantity (H) = E K +E V + E C + E D +E E is the HLF 
expectation of the Hamiltonian ifTOlfTTl with: 



Ek = J ?>)3^r^ x ' k )' 

Ey = J dx{7 tr (x)n(x), 
E C = 9 J dxn 2 (x), 



E D = - I dx$c(x)ri(x), 



rj f dxdk 

(2^)3 



$ B (x,k)W(x,k). 



(10) 

(11) 
(12) 
(13) 
(14) 



The contributions to (H) are the kinetic energy {En); the 
trap energy (£V); the combined direct and exchange con- 
tract interaction term {Ec)\ the direct dipolar term {Ed), with 
<&d (x) = J dx! t/dd(x — x')n(x'); the dipolar exchange in- 
teraction (Eg), where 



$ E (x,k) 



dk' 



f/ dd (k-k')W(x,k'), (15) 



with C/dd(k) = Cdd(cos 2 9u — |) the Fourier transform of 
t/dd( x ). The expressions for the interaction terms, ( fT2| i - ( |14) , 
are obtained using HF factorization to decompose second or- 
der correlation functions into products of single particle cor- 
relation functions, which can be expressed in terms of the 
Wigner function. Evaluating the above expressions within the 
HLF ansatz yields 



dx 



3k p (x) 3k z (x) 



U tT (x) + <?n(x) + -* D (x) + (x) 



K{x) (16a) 
1 n(x)l, (16b) 



with the local exchange term $ e (x) obtained from 

/(fir 
^$ B (x,k)^(x,k). (17) 

In addition to being local in position space, $e(x) has the 
simple analytic form 



$ £ (x) = C dd J[<5(x)]n(x), 



(18) 



where <5(x) = k z (x)/k p (x) — 1 is the relative distortion of 
the momentum distribution anc0 



J (it) = [yl + m (sinh 1 y/u)/y/u -l]/u- 



(19) 



is a monotonically decreasing function of u with J(0) = 0. 
Result ( [T8| shows that the effective exchange potential de- 
pends on the density and is only non-zero when the local 
momentum distribution is distorted from spherical symmetry 
[taking <5(x) = 0, $b(x) and Ee are zero and HLF reduces 
to Hartree theory]. The exchange potential appears with a pre- 
factor of T] in Eq. ( |16b| i and we find that S(x) > for bosons 
and <5(x) < for fermions so that Ee is always negative. 
The local form of exchange ( fT8] l we have arrived at is the cen- 
tral result that allows us to formulate a tractable and flexible 
theory. It is worth pausing to briefly compare to the HF treat- 
ment in which the full Wigner function needs to be evaluated 
and then convolved with the interaction potential to obtain the 
exchange potential ( p~5] > (e.g. see IfTOl QTJ). In contrast HLF 
theory does not require evaluating the Wigner function, yet 
contains the momentum dependence of the exchange term pa- 
rameterized by our two position dependent distortion fields 
{Kp(x), k z (x)} [or equivalently {m* (x), <5(x) }]. 

HLF equations: By requiring that SIhlf (|5]l is station- 
ary with respect to arbitrary variations of V c s(x), k p (x) and 
re z (x) we find: 



Kff(x) = f7 t x(x) + 2gn(x) + $ c (x) + ^(x), 
1 + §<5(x) 



rn 



(20) 
(21) 



m*(x) [l + 5(x)] 2 / 3 ' 

5(x) = j7? 4 dd [l+§f (x)][l+f (x)] J*[f(x)]. (22) 



! We note that (sinh 1 y^u)/y/u is real for u > — 1, that our J(u) 
+ u) 1/,3 ]/6 of 1211 . and that J(u) is easily differentiated for use 
(25} (see QU). 
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Equation ( 20 1 for the effective potential includes the local ex- 
change potential. The relative momentum distortion field <5(x) 
is determined by solving the transcendental Eq. p2|), and from 



2. 08 



this the effective mass is immediately given using Eq. (21 
We note that Eq. (21 1 ensures that the local kinetic energy is 



16a I is unity]. 



K(x) [i.e. the prefactor of K(n) in Eq. 

Equations (|20]>-([22|>, in conjunction with Eqs. Q and 
form the core set of equations of our theory that must be 
solved self-consistently. The direct potential, $d(x), can be 
efficiently computed using the convolution theorem. We note 
that a number of accurate and efficient techniques for doing 
this have been developed for the purpose of solving the Gross- 
Pitaevskii equation with DDIs (e.g. see [22]). 

HGF equations: We can develop a simplified version of 
HLF by setting a single global distortion, implemented by ig- 
noring the x dependence of the momentum distortion fields. 
Minimizing the free energy we find Eqs. (fT8j, ( |20] > and pi) 
(without position dependence of S or m*) and 



9?7Cdd / rfxn 2 (x) 
4~ J dxK{x) ' 



(1 + §£)(! + £)•/'(<*), (23) 



which we refer to as the Hartree Global-Fock (HGF) theory. 
The HGF theory captures the average exchange effects, and 
thus provides a good description of quantities such as the po- 
sition and momentum distributions. For many predictions the 
HGF theory will be inaccurate because the relevant properties 
are determined by local properties, e.g. mechanical stability 
is determined by the densest part of the gas near trap center, 
where local exchange effects are largest and drive the collapse 
to occur at lower dipole strengths. Similar considerations will 
be important in predicting phase transitions. HLF is just as 
easy to implement as HGF and calculation times are similar, 
with HGF approximately twice as fast as HLF calculations. 

Results: We validate the HLF theory by comparison to HF 
and Hartree calculations for a system in the harmonic trap 
t/ t r(x) = \mu: 2 (x 2 + y 2 + A 2 z 2 ), with A = uj z /uj p . To 
simplify our presentation we only discuss HGF calculations 
in cases that help illuminate its differences from HLF. 

For given /i and T we find that the HLF free energy is 
above, but close to the full HF value, and appreciably lower 
than the Hartree value. In Fig. [T] we compare the kinetic and 
dipolar-exchange energy (both give important contributions to 
il) for Bose and Fermi systems with fixed mean number of 
particles N. For the kinetic energy we find that HF and HLF 
calculations are in excellent agreement, and discernibly dif- 
ferent to the Hartree results. This difference, which is both 
positive and negative, arises directly from the momentum dis- 
tortion as well as from the self-consistent effects of interac- 
tions changing the chemical potential]^] 

The exchange energy Ee is zero for the Hartree theory and 
Fig.[TJb) shows the predictions of HF and HLF theories, again 
revealing excellent agreement. In the T — > limit of the 




FIG. 1. (color online) Comparison of energies for pure (g — 0) 
dipolar theories showing HF (crosses), HLF (solid curves), Hartree 
(dotted curves, zero for Ee so not shown), for aspect ratio A = 0.1 
(blue, dark grey), 1 (green, light grey), 10 (red, grey), (a) Ki- 
netic energy Ek, (b) dipolar exchange energy Ee- Main figures 
are for fermions with D t = 1 and insets are for a Bose gas at 
T = 1.5T,? . Dipole strength is parameterized in terms of Dt = 
CddiV 1//6 /(47rftu;aj5 ), where aho = y/h/mui and oj = ^/uipUi z , 

with T° = ^/N/C(3)Hw/k B and T% = \/Whu/k B the ideal 
gas Bose-Einstein condensation and Fermi temperature, respectively. 
Results in the insets terminate at finite D t due to approaching insta- 
bility (23). 

Fermi ga^jthe direct and exchange contributions are of sim- 
ilar magnitude (i.e. \Ee\ ~ \Ed\) for spherically symmetric 
traps ( A = 1), and \Ee\ is about an order of magnitude smaller 
than \Ed\ for the anisotropic cases with A = 10 and 0.1. 
This is because when the trap distorts the spatial distribution 



4 For T — > fermions in HLF we use n 
andK 



75 
:W 2 



{m*\ii- V cff ]} 3 / 2 /fi 3 , 



^|(m*) 3 / 2 [/x - V eB ] 5/2 /h 3 . This limit is difficult to realize 



' At fixed [i and T, the Hartree value for is less than HF and HLF. 



in HF calculations where the sharp Fermi surface in W(x, k) (e.g. see 
(9)) is smeared by the numerical grid resolution revealing an additional 
advantage of HLF. 
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away from being nearly spherical the direct interaction ( 13 1 is 
strongly enhanced, while \Ee\ remains roughly the same size. 
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FIG. 2. (color online) Comparison of pure (g = 0) dipolar theories 
showing HF (crosses), HLF (solid curves), HGF (dashed lines) and 
Hartree (dotted curves) in the radial (blue, dark grey) and axial (red, 
grey) directions, (a) <5(x) (Hartree results are not shown, being all 
zero) with A = 1 and Dt = 1 for bosons at T = l.2T° (upper 
curves) and fermions at T = O.lTp (lower curves), (b) n(x) and (c) 
momentum density n(k) = J W(x, k) for fermions at T = 

0.02T£ with A = 10 and D t = 2. The HGF agreement in (b) and 
(c) is reasonable and is not shown. 



The field <5(x), which is a key element of the HLF theory, 
is a measure of a local quadrupolar moment (i.e. distortion) of 
the momentum distribution, given by 



«5(x) 



IK (x) 
lk z (x) 



1, 7^(x) 



dk 



^ 2 VF(x,k), (24) 



where {7^ (x), 7/^ (x)} are the local momentum moments. 
We can use |24| to evaluate <5(x) from the full HF solutions 
[recall 5(x) — in the Hartree theory]. In Fig.[5Ja) we show 
the local momentum distortion for Bose and Fermi systems. 
Our results demonstrate that the momentum distortion varies 
spatially, with the largest distortion occurring at trap center 
(i.e. where density is highest), and that this effect is accurately 
captured by HLF theory. We also show HGF results which 
demonstrate that this approach predicts a reasonable average 
distortion, but fails to capture its spatial dependence. 

Both the Bose and Fermi systems exhibit similar behavior 
in their position-space distortion effects, i.e. the density elon- 
gates along the polarization (z) direction to reduce Ed. The 
momentum space distortion lfT2l is distinctive: to reduce Ee 
the Fermi system elongates along the k z direction whereas the 
Bose system reduces its k z extent to instead expand in the ra- 
dial momentum plane. This behavior is also apparent in the 
short range correlations between particles (e.g. see discussion 
in lfT2l ) and should be verifiable in current experiments 11241 . 
Density profiles in position and momentum space in Fig. |2|b) 
and (c), respectively, show that while the Hartree position den- 
sity is in reasonable agreement with HF and HLF, the Hartree 
theory fails to capture the difference between the momentum 
density in radial and axial directions. 

Conclusion and outlook: In this paper we have introduced 
a variational ansatz that converts the HF theory of dipolar 
Bose and Fermi gases to a local density dependent theory, with 
negligible error compared to the full HF solutions. The re- 
sulting calculations are practical to undertake with a dramatic 
reduction in required computing resources compared to HF 
and should support this burgeoning field of dipolar quantum 
gases. This approach also provides insight into the manifesta- 
tion of exchange effects in the dipolar gas, such as local and 
global momentum distortion, that could be verified in current 
experiments. In future work we will extend this approach to 
superfluid Bose and Fermi gases. 
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